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We investigate the stability of a spatially homogeneous and isotropic non-singular cosmological 
model. We show that the complete set of independent perturbations (the electric part of the per- 
turbed Weyl tensor and the perturbed shear) are regular and well behaved functions which have no 
divergences, contrary to previous claims in the literature. 
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I. INTRODUCTION 



The existence of singularities appears to be a property 
inherent to most of the physically relevant solutions of 
Einstein equations, in particular to all known up-to-date 
black hole and conventional cosmological solutions (Q). 
In the case of black holes, to avoid the singularity some 
models have been proposed(@, Q [J ) These models 
nonetheless are not exact solutions of Einstein equations 
since there are no physical sources associated to them. 
Many attempts try to solve this problem by modifying 
general relativity ([ESQ)- More recently it has been 
shown that in the framework of standard general relativ- 
ity it is possible to find spherically symmetric singularity- 
free solutions of the Einstein field equations that describe 
a regular black hole. The source of these solutions are 
generated by suitable nonlinear vector field Lagrangians, 
which in the weak field approximation become the usual 
linear Maxwell theory Similarly in Cosmol- 

ogy many non-singular cosmological models with bounce 
were constructed where the energy conditions or the va- 
lidity of Einstein gravity were violated. Such models are 
based on a variety of distinct mechanisms, such as a cos- 
mological constant (0), non-minimal couplings (H3), 
nonlinear Lagrangians involving quadratics terms in the 
curvature (|l4|), modifications of the geom etric structure 
of space-time quantum gravitv([l8j). and nonequi- 

librium thermodynamics (|l6j|). among others to restrict 
ourselves to homogeneous and isotropic solutions. Fur- 
ther investigations on regular cosmological solutions can 
be found in (E3)- 

In a previous paper(^i|) we have investigated a cos- 
mological model with a source produced by a nonlinear 
generalization of electrodynamics and succeeded to ob- 
tain a regular cosmological model. The Lagrangian of 
our model is a function of the field invariants up to sec- 
ond order. This modification is expected to be relevant 
when the fields reach large values, as is the case in the 
primeval era of our universe. The model is in the frame- 
work of the Einstein field equations and the bounce is 
possible because the singularity theorems are crr " 

cumvented by the appearance of a negative pressure (al- 
though the energy density is positive definite). Recently 



some papers started a detailed investigation of the tran- 
sition from contraction to expansion in the bounce of sev- 
eral models (|20|). In particular, in Einstein general rel- 
ativity, models with stress-energy sources constituted by 
a collection of perfect fluids and Friedmann-Robertson- 
Walker like geometry were examined(|24|). The claim in 
that paper is that a generic result about the behavior of 
scalar adiabatic perturbations was obtained. The result 
is the following: scalar adiabatic perturbations can grow 
without limit in two situations represented by the points 
where the scale factor attains its minimum value and 
where p + p = 0. The first point corresponds to the mo- 
ment in which the Universe passes through the bounce; 
the second corresponds to the transition from the region 
where the Null Energy Condition (NEC) is violated to 
the region where it is not. We will show that these in- 
stabilities are not an intrinsic pro perty of a model with 
bounce as claimed in reference |24) but a consequence of 
the existence of a divergence already appearing in the 
background solution if the description of the source is 
made in terms of a perfect fluid. We will present a specific 
example of a model with bounce, generated by a source 
representing two non-interacting perfect fluids, that has 
regular perturbations in the situations described on ref- 
erence 1241. 



II. THE MODEL 

We limit our analysis to a model 0] in which the 
singularity of the Friedmann-Robertson- Walker (FRW) 
geometry is avoided by the introduction of nonlinear cor- 
rections to Maxwell electrodynamics. We will consider 
modifications up to second order terms in the field in- 
variants 



L 



1 



F + aF 2 + f3G 2 , 



(1) 



where F = F lw F' i \G = \r) aSitlv F ati F^ ', a and /3 are ar- 
bitrary constants|2q. The term FG will not be included 
in order to preserve parity. The energy-momentum ten- 
sor for nonlinear electrodynamic theories reads 



-4 Lp F^ a F au + (GLq — L) 9fiu , 



(2) 
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where Lp represents the partial derivative of the La- 
grangian with respect to the invariant F and similarly 
for the invariant G. In the early universe, matter should 
be identified with a primordial hot plasma |29l l3G| and 
as a consequence we are led to the case in which only the 
average of the squared magnetic field survives (see 0] 
and references cited there for details). Since the average 
procedure is independent of the equations of the electro- 
magnetic field we can use it in the generic expression of 
the energy-momentum tensor to obtain 



T MV = (p + p) v^Vv -pg^, 



where 



Pi 



Pi 



H 2 (1 ~8aH 2 ), 



H 2 (l-A0aH 2 



(3) 

(4) 
(5) 



The standard result of the linear Maxwell theory can be 
recovered by setting a = f3 = 0. 

We set for the fundamental line element 

ds 2 = di 1 - ^® [dr 2 + r 2 (d9 2 + sin 2 9 dip 2 )) , (6) 

where e = —1, 0, +1 hold for the open, flat (or Eu- 
clidean) and closed cases, respectively. The Einstein's 
equations and the equation of energy conservation writ- 
ten for this metric become: 



,o 2 e 1 



P'j + HPl +P~f)- 



(7) 



(8) 



(9) 



Inserting Q and © in yields for the magnetic field: 



H = Hna 



(10) 



where H Q is an arbitrary constant. With this result the 
equation J7J) can be integrated. For the case e = the 
solution is: 



a{tf = Ho 



\{kc 2 t 2 + 12a) 



1/2 



(11) 



The interpretation of the source as a one component 
perfect fluid in an adiabatic regime has some difficulties 
that are at the origin of the instabilities found in f[2^1. 
The sound velocity of the fluid in this case is given by 

(M) 



dp 1 



Pj_ 

Pi 



Pi 



0(Pi + Pi) 



(12) 



This expression, involving only the background, is not 
defined at the points where the energy density attains an 
extremum given by 9 = and p 7 +p 7 = 0. In terms of the 
cosmological time they are determined by t = and t = 
±t c = 12a/fcc 2 . These points are well-behaved regular 
points of the geometry indicating that the description of 
the source is note appropriate. 

This difficulty can be circumvented if we adopt an- 
other description for the source of the model. This can 
be achieved if we separate the part of the source re- 
lated to Maxwell dynamics from the additional non- linear 
a— dependent term on the Lagrangian. By doing this the 
source automatically splits in two noninteracting perfect 
fluids: 



T, 



where, 



T lv = (Pi + Pi) v H v v - Pi 9nu, 



T lu = (P2 + Pi) V^V V - p 2 SV, 



and 



Pi 
Pi 

P2 
P2 



20 „ 4 
--aH 



(13) 

(14) 
(15) 

(16) 

(17) 
(18) 
(19) 



Using the above decomposition it follows that each one 
of the two components of the fluid satisfy independently 
equation This indicates that the source can be de- 
scribed by two non-interacting perfect fluids with equa- 
tion of states pi = 1/3 pi and P2 — 5/3/92- The equation 
of state for the second fluid should be understood only 
formally as a mathematical device to allow for a fluid 
description. 



III. GAUGE INVARIANT TREATMENT OF 
PERTURBATION 

In a series of papers 0, |3^, |3i| we have established 
a complete and self-consistent theory to deal with the 
problem of perturbations of the FRW cosmology. The 
very well-known problem of the gauge dependence of the 
perturbations was addressed and solved by the introduc- 
tion of a complete set of gauge invariants variables that 
represents direct observable quantities. We present here 
a summary of this formalism in order to fix the notation 
and to aim for the self consistence of this paper. 

The source of the background geometry is represented 
by two-fluids, each one having an independent equation 



3 



of state relating the pressure and the energy density 
(j>i = Xipi, where the i = 1,2. Following the standard 
procedure we consider arbitrary perturbations that pre- 
serves each equation of state. Thus the general form of 
the perturbed energy-momentum tensor is written as 



(1 + Aj) S(piVfj,v u ) — XiSiflig^)- (20) 



The background geometry is conformally flat. Thus 
any perturbation of the Weyl tensor is a true perturba- 
tion of the gravitational field. It is convenient to repre- 
sent the Weyl tensor W a pys in terms of its corresponding- 
electric and magnetic parts (these names come from the 
analogy with the electromagnetic field) : 



E aP = -W a ^ v v"v v 



H, 



a/3 



-W, 



(21) 



(22) 



These variables have the advantage that since they are 
null in the background, their perturbations are gauge in- 
variant quantities j34[. These definitions imply for the 
tensors and the following properties: 



E 



E vll , E„ v v» = 0, E^ v = 0, 



and for the magnetic tensor: 



0, H^g 



0. 



A. Some mathematical machinery 

The metric and the vector v" (tangent to a time- 
like congruence of curves T) induce a projector tensor 
h^ v which separates any tensor in terms of quantities 
perpendicular and parallel to v". The projector is defined 



7 A r A A 



and has the property that 



h h vX — h ; 



(23) 
(24) 



The equations of motion for the first order perturba- 
tions are linear so it is useful to develop all perturbed 
quantities in the spherical harmonics basis. In this paper 
we will limit our analysis to perturbations represented in 
the scalar base defined by the equation: 



2 1 n 



(25) 



where m is the wave number and \7u is the covariant 
derivative in the hypersurface with normal v" and metric 
h^v. From now on we will suppress the index n. 

The scalar Q allows us to define the associated vector 
7Tu and traceless tensor Pu V '■ 



(26) 



(27) 



In the case of scalar perturbations the fundamental set 
of equations, determining the dynamics of the perturba- 
tions are (see the appendix): 



(SE^yh^hu + {5E» v yh^hf + 0(5E° P + 5E° 2 



? af3 



a0 



(SE 1 ^ + SEl^h^ = -(6 Pl + S P2 ), a h° 
- gPi 5vf - -p 2 5v e 2 



(28) 



(29) 



(6al. f) + dal. f} ) h^ a Kf 



+ 3© [Saj^ + 5*%) = -6 El - SE^O) 



-XxS(p,ph p ll ) + (l + Xi)p5al = (31) 

- X 2 S ( PJJ h r \) + (1 + X 2 )p Sal = 0. (32) 

The acceleration a", the expansion O and the shear a^ v 
that appear in the above equations are parts of the ir- 
reducible components of the covariant derivative of the 
velocity field defined as: 



a 1 " = v^.,,v 



V U 1 



e = «* 



= -h" {a h )[j v v i 



(33) 
(34) 

(35) 



The expansion of the pert urbations in terms of the spher- 
ical harmonic basis is 



: perti 

s 27] 

6p = N(t)Q, 



Sv^ = V(t)V ia Q la , 
5a" = VU ia Q\ a: 
SE" V = E{t)P^ , 
8<r»" = T.(t)P^ v , 



(36) 
(37) 
(38) 
(39) 
(40) 
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B. Perturbation of a bouncing universe 

After presenting the necessary formalism we shall to 
start the analysis of the perturbations of the bouncing 
cosmological model displayed in previous section. Using 
the above expansion into the equations (cf. apendix) 
(fTUOf t. OH}, GJE), EH) and (Eg we obtain: 



Ei + i?2 



6e + k- 



(Ni + N 2 - piVi - P2V2) , (41) 



We should analyze the behavior of these perturbations 
in the neighborhood of the points where the energy den- 
sity attain an extremum. This means not only the bounc- 
ing point but also the point in which p + p vanishes. Let 
us start by the exam at the bouncing point t = 0. 

The expansion of the equation of Ex in the neighbor- 
hood of the bouncing, up to second order, is given by: 



Ei + actEx + (b + bit 2 )E 1 = 



(52) 



The constant a and the parameter b are defined as follows 



Ei+E 2 + -e(Ex + E 2 ) 



1 + Xx\ ( 1 + Ai . 

- ( — - — ] p 2 £ 2 (42) 



Si 



Vx — V2 = —Ex — E 2 , 



(43) 



2fi 



V6H t c 



(53) 



(54) 



Ai (Nx-pxVx) + (l + Xx)piVx = 0, 



(44) 



, = __6 3_ 

1 2tl 4tf. 



(55) 



A2 (iVjj — P2V2) + (1 + A2) P2V2 = 0, (45) Defining a new variable f as 



These equations can be rewritten in a more convenient 
way as: 

/2Ai(3e + fc 2 ) 



Si 



Va 2 (l + Ai)pi 



1 £1 



(46) 



/(*) = ^(tjexp +^ - -sjbx - ^ t 2 , (56) 



doing a coordinate transformation for time as indicated 
bellow 



Z = -i(\lbx- ^)t\ 



(57) 



Ex + ^QEx = -\{l + \x)pxLVx, (48) 



£ 2 + ie£ 2 = -i(l + A 2 ) /02 E 2 , (49) 

As we have shown in (|3l|), the whole set of scalar per- 
turbations can be expressed in terms of the two basic 
variables: Ei and The corresponding equations can 
be decoupled. The result in terms of variables Ei is the 
following: 

h + i±^ + { ^e 2 



We obtain the following confluent hypergeometric equa- 
tion (H|) 



where 



e/+(l/2-0/ + e/ = 0, 
i(b-a/2) 1 



' 4(6i -a 2 /4)V2 2' 
The solution of this equation is given by: 



f(t) = A M(d, 1/2, -i(\ bx--r)t') 



(58) 



(59) 



(60) 



(3 + 



(50) 



Note that there is no sum in indices and j 7^ i in this 
expression. In our case the values of Aj are Aj = (|, |). 
In the first case the equation for the variable Ei became: 



Ei + -QEx 



l n2 5fc 2 
3 9 - pl - p2 ~M 



E {1) = (51) 



where A is an arbitrary constant and M(d, 1/2, £) con- 
fluent hypergeometric function. The confluent hypergeo- 
metric function is well behaved in this neighborhood and 
so also is the perturbation i?i(f) given by: 



Ex(t) = Re[AM(d,l/2,-i(\/bx-—)t') 



* exp _ 4 + 2V 6i "Tr ] - 



(61) 
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The perturbation E 2 at this same neighborhood, after 
the same procedure we did before result in the following 
equation: 



E 2 + atE 2 + {b + b 1 t 2 )E 2 = 



(62) 



This is the same equation we obtained for E\ , they differ 
only by the values of the parameters a, b and b± that in 
this case are: 



9 

2i! 



2/? 



V6H t c 



bi = 



5m 2 



t 3 c H Ve ti 

Then the solution in this case is: 



(63) 



(64) 



(65) 



E 2 (t) = Re[AM[ d,l/2,-i{\jb x - — \t 2 ) 



* ex p-\i~2\l bi ~T) t2] 



(66) 



Again the confluent hypergeometric function is well be- 
haved in this neighborhood and so also is the perturba- 
tion E 2 (t). At the neighborhood of the point t = t c the 
equation for the perturbation E\ is given by 

E x + aEi + {b + bit) E 1= (67) 

where the parameters a, 6, 61 in this case are given by 

5 



a = 



4f, 



3 V3m 2 



V3 ( m 2 3 



At 2 \3H 2t c 



(68) 



(69) 



(70) 



We would like to remark that this equation is different 
from the equations (|52I62|) obtained in the neighborhood 
of t = 0. We proceed doing the following variable trans- 
formation: 



E x (t) = cxp~—w(t) 



(71) 



The differential equation for this new variable is 

w+(b-{a/2) 2 + b 1 t)w = (72) 



The solution for this equation is 

b- {a/2) 2 



w(t) 



wq AiryAi 



bit 



b 2/3 



(73) 



The AiryAi are regular well behaved functions in this 
neighborhood and so also the perturbations E x . Finally 
we look for the equation of E 2 at the neighborhood of 
t = to, it becomes 

E 2 + aE 2 + (b + ht) E 2 = (74) 

where the parameters a, 6, b x in this case are given by 

(75) 



9 

4^ 



b = 



t c \ At c 



V3m 2 
6H 



h = 



5V3 ( 1 



2t\ \t c 6H a 



(76) 



(77) 



This equation differ from ea. (|67|) only by the numerical 
values of the parameters a, b, b x so we obtain the same 
regular solution 



/ at A , r b- (a/2) 2 +M, 
E 2 = Re I exp -— w Q AiryAi[ ] 



IV. CONCLUSION 



(78) 



Recently there has been a renewed interest on nonsin- 
gular cosmology. As a direct consequence of this some 
authors have argued a gain st these models based on in- 
stability reasons. In ( 24]) it has been argued that a 
rather general analysis shows that there are instabilities 
associated to some special points of the geometrical con- 
figuration. They correspond to the points of bouncing of 
the model and maxima of the energy density, where the 
description of the matter content in terms of a single per- 
fect fluid does not apply. In the present paper we have 
shown, by a direct analysis of a specific nonsingular uni- 
verse, that the result claimed in the quoted paper does 
not apply to our model. We took the example from a re- 
cent paper ( [Tgj | ) in which the avoidance of the singularity 
comes from a non linear electrodynamic theory. We used 
the quasi Maxwellian equations of motion ([3l|], |23, |33|) 
in order to undertake the analysis of the perturbed set 
of Einstein equations of motion. We showed that in the 
neighborhood of the special points in which a change of 
regime occurs, all independent perturbed quantities are 
well behaved. Consequently the model does not present 
any difficulty concerning its instability. This paves the 
way to investigate models with bounce in more detail 
and to consider them as good candidates to describe the 
evolution of the Universe. 
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V. APPENDIX A: COMPARISON WITH 
OTHERS GAUGE-INVARIANT VARIABLES 

FRW cosmology is characterized by the homogeneity 
of the fundamental variables that specify its kinematics 
(the expansion factor 6), its dynamics (the energy den- 
sity p) and its associated geometry (the scalar of curva- 
ture R). This means that these three quantities depend 
only on the global time t, characterized by the hypersur- 
faces of homogeneity. We can thus use this fact to define 
in a trivial way 3-tensor associated quantities, which van- 
ish in this geometry, and look for its corresponding non- 
identically vanishing perturbation. The simplest way to 
do this is just to let U be a homogeneous variable (in the 
present case, it can be any one of the quantities p, or 
R), that is U = U(t). Then use the 3-gradient operator 
(3) Vm defined by 



factor 0. This allows us to define the quantity r\ a as: 



(3) 



Va 



to produce the desired associated variable 



(79) 



(80) 



In [3|j these quantities were discussed and its associ- 
ated evolution analysed. In the present section we will 
exhibit the relation of these variables to our fundamen- 
tal ones. We shall see that under the conditions of our 
analysis |39j these quantities are functionals of our basic 
variables (E and E) and of the background ones. 



(84) 



Using the constraint relation eq. l|97|l we can relate this 
quantity to the basic ones: 

= (l-^) 0i . (85) 



C. The geometrical variable r 

We can choose the scalar of curvature R which depends 
only on the cosmical time t like p and to be the U- 
geometrical variable. However it seems more appealing 
to use a combined expression r involving R, p and 
given by 



t = R + (1 + 3A) p 



0- 



(86) 



In the unperturbed FRW background this quantity is de- 
fined in terms of the curvature scalar of the 3-dimcnsional 
space and the scale factor A(t): 

&R 
A 2 ' 

We define then the new associated variable t„ as 



A. The matter variable Xi 

It seems useful to define the fractional gradient of the 
energy density Xa as |3(J 



1 



(3) 



V« P- 



(81) 



Such quantity x.a is nothing but a combination of the 
acceleration and the divergence of the anisotropic stress. 
Indeed, from the above equations it follows (in the frame 
in which there is no heat flux) 



(87) 



This quantity r Q vanishes in the background. Its pertur- 
bation can be written in terms of the previous variations, 
since Einstein's equations give 



1 



VI. 



0' 



APPENDIX B: QUASI-MAXWELLIAN 
EQUATIONS 



S X i = Sat + Y p ^ 



(82) 



/From what we have learned above it follows that this 
quantity can be reduced to a functional of the basic quan- 
tities of perturbation, that is S and E, yielding 

^■- 2 ( 1 -?)^( E -| E )«- (83 » 



B. The kinematical variable rji 

The only non-vanishing quantity of the kinematics of 
the cosmic background fluid is the (Hubble) expansion 



We list below the quasi-Maxwellian equations of grav- 
ity. They are obtained from Bianchi identities as true 
dynamical equations which describe the propagation of 
gravitational disturbances. Making use of Einstein's 
equations and the definition of Weyl tensor, Bianchi iden- 
tities can be written in an equivalent form as 

2 6 y 



Using the decomposition of Weyl tensor in terms of 
E a p and H a p (see Section IHIjl and projecting appro- 
priately, Einstein's equations can be written in a form 
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which is similar to Maxwell's equations. There arc 4 in- 
dependent projections for the divergence of Weyl tensor, 
namely: 

W a ^. v VpV„ h a \ 
W a ^ ;v rf x a p V^V X , 
W a ^ ;v ^ x a(j V X , 

The unperturbed quasi-Maxwellian equations are thus 
given by: 

h ea h x ~< E aX - n + V s ^ V V H vX a» x + 3H £V u v 

+ \*?» <V + \h ea ^"V (88) 

h ea h x t H aX;7 - rf ppa , V? E vX a^ x - 3E™ u u 
= (P + PK -\r! ea(ix Vxq a ;p 

+ |f^(^+^)/ ^ (89) 



The contracted Bianchi identities and Einstein's equa- 
tions give the conservation law 

T» v . v = 0. 

Projecting it both in the parallel and the orthogonal sub- 
spaces we obtain: 

T^,,V a = 0, 

V = o, 

which give the following equations: 

p + ( P + p)e + n + q a . a - = o, (92) 

(p + p)a a - p. a h» a + q^ a + Oq a 

+ q v Q av + q v u av + it a v . u 

+ 7r^Q^V a = 0, (93) 

and, from the definition of Riemann curvature tensor 

we obtain the equations of motion for the unperturbed 
kinematical quantities as: 

e + ^ + + 2uj 2 -a a ;a = R^V»V, (94) 



h„ s h v x + <dH eX - ^H^ {£ h x ^ V^ v 

+ r^^rff*" VfjVr H ai Q v p 

- a a E^rf* af> y 7 

4. I ErS h ( e T 7 A)7Q ' 3 V 

= -\q (e u x) + \h eX tfuip 
+ \<?p {e v X)a ^v a q a 

+ \h v ^ a ^ V„ n va;P (90) 

h» e h x E^ + OE eX - ^E u {e h x \ V*" 
+ r) Xv ^rfV ra V a V T E ai Q vP 

- \(p + p> eX + \q ( ^ X) 

- \h^h x ^ q Ka + l -K e h a x ^ 

+ ^(V^-iTr^^ + ie^l) 



W ^ + ^h a p(-2uj A - 2a A + 

a X -x)+a a a [} - -h^hp" (a^ v + a„ ;fl ) 
2 

+ ^^Q/J + ^au^p + ^au^p 

= RaeflvV'V - \R^V V Kp,{Qh) 

1 2 

h^hp" Qftv - -h a tl h f3 v {a ll . v - a v , u ) + -9u> Q/3 

+ CT a ^p - ap n UJ^a = 0. (96) 

We also obtain from the definition of R a p U v three con- 
straint equations: 

^Q^x-i^i+^-iUWx-a^axv+uxv) = R^V»h» x , 

(97) 

(98) 



w% + 2uj a a a = 0, 



- 1 K e h x a f]/ 11 ' V u (<j a p + uj a p)- n + «(r U\) = H tX . 

(99) 

These results constitute a set of 12 equations which 
will be used to describe the evolution of small pertur- 
bations in FRW background. Writing all the perturbed 
quantities in the form 

X (perturbed) (background) ~t~ 
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and after straightforward manipulations we finally obtain (5H, V ) = — — h a h 13 )((<5cr Q7 ) \ + (5u) a y) a) f]f) eiX Ve 
the perturbed equations from the set of equations JSSJl- 2 

m as: ( 109 ) 

„ 1 , « (Sp)' + Q(S P + Sp) + { P + p) (60) + (Sq a )., a = (110) 

+ f ^ E ^ Q7rA W(^ eA ) V P (<%) + P,o (^°) - (5p) >j8 ^ + (p + p) (6a,) 

_ ItfHx"). h ( a v ^ x V T + V«(*0* + 3 (%) + V» ( 5 ^ af> ):,0 = 0. 

2 1 "' ' " (HI) 

= _-(p + p)(^) 

+ g ^ - \ h* a hP" (%);, 

+ i h^ a h > (su^y + 1 e (<m Q/3 ) 
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(iff"*')* ft^/i/ + 9 {6H af3 ) - ]^{8H v (a )h (3 \ V*" 
+ f V M£ r, aXT1 V M V T (SH ej ) h Xv 

= \ V<V )eT " V, (m vs ). T (101) 



(6H atl ). v h as W = (p + P ) (6u> s ) - l - rf<*^ V, (5q a ), 

(102) 

[8E ail ), v h as h^ v = ^(Sp), a h as ~^p(SV") 



(sey + (se) - (Sa a \ a = _Ii±^) (d>) (104) 

(Scr^y + -/i M „((5a Q ) ;Q - -(£a( Q ) ; /3) ^" h ^ 



r ^6 (5(7^) = -(6E^) - (105) 



+ |e - ir, Q ^ (« a/J ) n V Q (106) 



jj(*e), A fc\, - ~e (sv^) + 1 e (sv°) 

- (^ + ^). a ^ = -(*?„](107) 



(<W a ) ;ct = (108) 
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